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FOREWORD 


This  technical  report  is  written  in  the  course  of  a program  entitled: 
"Linearized  Theory  of  Gas-Dynamic-Laser  Nozzle  Wakes  and  Applications," 
sponsored  by  AFSWC,  Kirtland  AFB,  New  Mexico  and  monitored  by  AFOSR  under 
Contract  No.  F44620-75-C-0016.  This  phase  of  the  work  deals  with  the  general 
subject  of  the  flow  exiting  adjacent  supersonic  nozzles.  Specifically,  it  is 
sought  to  characterize  and  predict  the  development  of  wakes  emanating  from 
the  trailing  edges  of  the  cusps  formed  by  such  adjacent  nozzles.  Wake  meas- 
urements are  carried  out  at  the  Aeronutronic  Ford  Supersonic  Wind  Tunnel 
under  this  program;  the  data  will  then  be  used  to  test  the  adequacy  of 
theories  predicting  the  wake  development.  This  report  summarizes  the  basic, 
linearized  flow  theory  predicting  the  development  of  the  wake,  for  comparison 
with  the  experimental  data. 

The  author  is  indebted  to  P.J.  Ortwerth  of  AFSWC  for  guidance  and  encourage- 
ment in  this  work,  and  to  Milton  Rogers  of  AFOSR  for  his  assistance  in  this 
program. 
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ABSTRACT 


The  objective  of  this  report  is  to  assemble,  in  compact  form,  the  equations 
describing  the  development  of  wakes  trailing  from  gas-dynamic-laser  nozzle 
cusps.  A particular  family  of  appropriately  short,  geometrically  similar 
DeLaval  nozzles  are  considered.  The  exit  flow  is  assumed  parallel  and  the 
trailing  edge  thickness  is  taken  to  be  sharp.  Fully  laminar  and  fully 
turbulent  wakes  are  considered  for  arbitrary  exit  Mach  and  Reynolds  number, 
specific  heat  ratio  and  nozzle-to-stagnation  temperature  ratio.  The  total 
viscous  drag  and  heat  exchanged  with  the  nozzle  cusp  are  calculated,  for 
such  conditions,  for  both  laminar  and  turbulent  nozzle  boundary  layer.  To 
calculate  the  laminar  wake,  Kubota’s  linearized  similarity  theory  is  then 
utilized;  for  the  turbulent  wake  an  analogous  similarity  theory  is  used, 
with  key  information  on  turbulent  transport  drawn  from  recent  two-dimensional 
supersonic  wake  experiments.  The  results  are  given  in  non-dimensional 
form  by  expressing  distances  in  terms  of  the  nozzle  exit  height.  Most  wake 
parameters  are  found  insensitive  to  Reynolds  number  but  scrongly  dependent 
on  the  other  design  parameters.  Procedures  are  indicated  and  charts  and 
tables  are  supplied  by  which  the  reader  can  apply  the  formulas  to  any  set 
of  design  parameters,  and  a specific  example  at  Mach  4 is  computed  in  detail. 
Assumptions  and  approximations  are  pointed  out  for  future  test  by  experiment. 
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I . INTRODUCTION 


The  technique  of  "stacking"  supersonic  DeLavai  nozzles  on  top  of  one  another 
is  at  present  the  standard  method  of  expanding  a flowing  gas  to  low  temper- 
ature supersonic  flow  over  a short  physical  distance.  This  "nozzle  array" 
allows,  in  principle,  the  expansion  length  to  be  small  regardless  of  the 
mass  flow,  and  is  therefore  a good  scheme  for  high-power  gas-dynamic  lasers. 
Ideally,  uniform  supersonic  flow  is  obtained  in  the  discharge  volume  (the 
"laser  cavity").  In  practice,  however,  each  nozzle  cusp  sheds  a supersonic 
wake  which,  by  definition,  introduces  flow  non-uniformity  in  the  cavity. 
Non-uniformities  degrade  the  laser  performance,  and  it  is  thus  important  to 
predict  the  cusp  wake  flow  for  each  given  set  of  GDL  operating  conditions. 

This  report  presents  formulas  describing  laminar  and  turbulent  wake  behavior 
for  various  nozzle  exit  Mach  and  Reynolds  numbers,  specific  heat  ratios  and 
nozzle  wall  temperatures.  Such  generality  is  achieved  by  'Standardizing" 
the  nozzle  contour  and  utilizing  the  linearized,  similarity  approach  to  wake 
flow.  For  laminar  wakes  this  linearized  solution  has  been  available  for  some 
time.  A corresponding  solution  has  been  available  for  the  turbulent  wake, 
with  the  transport  coefficients  given  from  experiments  on  two-dimensional 
supersonic  wakes.  The  disadvantage  of  the  linearized  solution  is  its  inval- 
idity near  the  cusp  trailing  edge,  and  also  when  the  design  parameters  lie 
on  certain  extreme  ranges.  The  advantage  lies  in  the  compactness  with  which 
the  results  are  expressed  in  charts  and  tables.  Also,  the  qualitative  trends 
in  wake  behavior  with  changing  parameter  values  should  not  then  change  much 
if  non-linear,  exact  solutions  were  to  be  found. 

2.  STATEMENT  OF  THE  PROBLEM  AND  SUMMARY  OF  THE  SOLUTION 

The  geometrical  configuration  of  the  problem  is  shown  on  Fig.  1.  A nozzle 
"cusp"  is  enclosed  symmetrically  by  two  nozzle  contours.  The  nozzle  exit 
flow  is  parallel  to  the  axis  and  the  cusp  T.E.  (trailing  edge)  is  taken  as 
infinitely  thin.  The  nozzles  are  assumed  perfectly  two-dimensional,  i.e., 
no  interaction  is  allowed  between  their  flow  and  the  boundary  layer  growing 
on  the  nozzle  sidewall.  Two-dimensionality  is  also  assumed  for  the  wake  and 
nozzle  flow  downstream  of  the  T.E. 

Because  of  momentum  and  heat  balance,  the  wake  integral  properties  are  con- 
trolled by  what  happens  on  the  nozzle  surfaces.  Therefore,  it  is  first  im- 
portant to  "standardize"  the  nozzle  contour,  i.e.,  to  find  a family  of  nozzles 
whose  non-dimensional  coordinates  are  functions  of  M^  alone;  this  is  done 
in  Section  5.  Then,  Section  6 calculates  the  momentum  and  heat  exchange 
between  nozzles  and  flow,  as  a function  of  M^Re',  y and  TW/TQ.  This  cal- 
culation is  done  both  for  laminar  and  turbulent  flow  (Sections  7.1  and  7.2 
respectively. ) 

In  Section  7.1,  the  linearized  laminar  wake  theory  of  Kubota's  is  written 
down  and,  with  the  inputs  from  Section  5,  used  to  express  specific  variations 
for  the  laminar  wake  properties.  The  same  is  done  for  turbulent  wakes  in 
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Section  7.2,  except  that  Kubota's  theory  is  replaced  by  the  standard  linear- 
ized turbulent  wake  analysis  with  pivotal  inputs  on  turbulent  transport 
coefficients  taken  from  experiment. 


Sections  7.1  and  7.2  contain  all  the  formulas  needed  to  predict  the  longi- 
tudinal and  lateral  variations  of  fluid  velocity,  temperature,  etc.,  as  a 
function  of  the  GDL  design  parameters.  A guide  for  using  these  formulas  in 
step-by-step  fashion  is  presented  in  Section  8.1,  and  a numerical  example  of 
this  step-by-step  process  is  given  in  Section  8.2.  Finally,  a complete  set 
of  computations  is  given  in  Section  8.3  for  the  M,^  = 4,  y = 1.4  case,  which 
is  the  case  selected  for  experimental  verification  in  the  Aeronutronic  Ford 
Supersonic  Wind  Tunnel. 

3.  ASSUMPTIONS  AND  APPROXIMATIONS 

The  following  assumptions  and  approximations  will  be  made: 

A.  The  subsonic  portion  of  the  nozzle  cusps  contributes  nothing  to  the 
total  nozzle  heat  transfer  and  skin  friction. 

B.  The  friction  and  heat  transfer  coefficients  at  any  distance  p (measured 
from  the  throat  along  the  nozzle  axis)on  the  nozzle  surface  are  equal 
to  those  at  distance  p from  the  leading  edge  of  a flat  plate  in  uniform 
flow,  at  edge  conditions  equal  to  those  on  the  nozzle  axis  at  £,  and 
without  pressure  gradient. 

C.  The  laminar  friction  and  heat  transfer  coefficients  are  independent 
of  the  local  Mach  M,  and  of  the  local  heat  transfer  rate. 

D.  The  turbulent  friction  coefficient  is  independent  of  heat  transfer. 

E.  The  Prandtl  number  along  the  nozzle  wall  is  constant. 

Note  that  assumption  (A)  is  especially  reasonable  for  GDL  nozzles,  which  are 
almost  always  designed  to  contract  and  expand  the  flow  quickly.  Fig.  1 shows 
typical  GDL  nozzle  geometries. 

4.  GENERAL  EXPRESSIONS  FOR  GDL  NOZZLE  FRICTION  FORCE  AND  HEAT  EXCHANGE 

In  Sections  4-6,  we  use  the  coordinate  £ to  describe  streamwise  distance 
along  the  nozzle  axis,  with  £ = 0 at  the  nozzle  throat  (see  Fig.  1).  In 
later  Sections,  the  coordinate  x,  zero  at  the  T.E.,  is  used  in  connection 
with  the  wake  flow. 

Since  the  flow  is  two-dimensional,  the  total  friction  force  per  unit  span 
on  the  nozzle  "cusp"  (Fig.  1)  is  twice  the  integral  of  the  shear  force  from 
the  throat  (|  = 0)  to  the  trailing  edge  (T.E.): 
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J 


£ = 7 

Y CO  J Co 


u u / 


Since  M is  a function  of  e/i,  in  general,  our  objective  is  in  principle  at- 
tainable since,  upon  integration,  the  integral  in  eq.  (5)  or  eq.  (6)  will 
only  be  a function  of  exit  quantities  (subscript  «) . 

5.  PRESCRIPTION  OF  THE  GDL  NOZZLE  CONTOUR 

As  just  mentioned,  the  integral  in  eq.  (5)  and  eq.  (6)  can  be  integrated  into 
an  expression  involving  nozzle  exit  quantities  only.  For  any  given  nozzle, 
the  function  M = M(^/^,)  is  also  given,  by  definition.  Then,  the  quantities 
in  the  integrand  can  be  expressed  in  terms  of  M and  thus  also  in  terms  of 
^/i;  for  example: 


£ = Z (Miln  = 1.1) 
Pc  p‘  PoV 


This  procedure  will  be  displayed  in  detail  later,  when  the  integrals  (5) 
and  (6)  are  actually  computed  (see  Section  6.) 

The  above  procedure,  however,  cannot  be  generalized  for  arbitrary  nozzles  of 
arbitrary  exit  M^,,  since  it  is  well  known  that  there  is  no  unique  function 
M(?/<t)  for  DeLaval  nozzles.  We  wish  to  set  as  an  objective  the  evaluation 
of  the  integrals  (5^  and  (6)  for  such  "general"  nozzles  and  for  arbitrary 
Mm.  Thus  we  week  a function,  say  z(M),  which  depends  on  £/£  so  that  does 
not  appear  explicitly.  One  of  the  simplest  choices  is 

? = 2(M)  = ^ = 0 M = l ( 5/H  = 0)  (12) 


Mof  * 


= I @ M=  Mqq  l)  (13) 


iqaawMw— iui  i ijumw  >w  warn 


Consider  now,  that  there  is  a function  z(^/^,)  which  is  unique  for  all  DeLaval 
nozzles.  Then  for  each  a z is  obtained  and,  through  eq.  (14),  also  a M 
once  the  exit  is  given.  The  strategy  of  evaluating  the  integral  in  (5) 
or  (6)  once  M,*,  is  given  is:  first,  express  all  quantities  in  the  integrand 
in  terms  of  M,  then  in  terms  of  z (using  eq.  (12)  ) and  then  in  terms  of 
using  the  known  z(^/<t).  The  integral  can  then  be  evaluated,  the  final 
expression  containing  (and  also  Re«,  Y,  etc.,  as  they  appeared  in  the 
integrand).  Thus  the  viscous  force  of  the  nozzle  cusp 


-—-■L  = / ^eoo  » Y *fc.)  (15) 

2.  PcDU® 


As  mentioned  above,  there  appears  to  be  no  unique  function  z(|/^)  for  all 
nozzles.  However,  Shames  and  Seashore  (Ref.  1)  indicate  that  there  exists 
a family  of  "short"  supersonic  nozzles,  treated  in  their  paper,  which  shows 
distinct  geometric  similarity  for  all  M^.  Such  "short"  nozzles  are  obviously 
suitable  for  GDL's.  Availing  ourselves  of  this  coincidence,  we  make  the 
following  assumption  in  this  work: 

(F)  "The  supersonic  nozzles  discussed  here  (also  called  "GDL" 
nozzles)  are  characterized  by  a unique  function 
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(16) 


The  next  step  is  to  find  z(^/-f-).  A GDL  nozzle  for  = 4 has  been  designed, 
and  its  M at  each  £ (or  ^/-t)  is  shown  on  Table  I.  From  M and  the  variation 
of  z vs  £/£  was  plotted  on  Fig.  2.  This  curve  was  next  least-squares  fitted 
by 


2 = 


(17) 


Trials  were  made  with  polynomials  of  degrees  4 through  8,  the  coefficients 
being  shown  on  Table  II.  The  8th  degree  polynomial  best  fits  the  z(^/-t) 
curve,  as  shown  on  Fig.  2,  but  all  integrals  were  also  evaluated  using 
lesser-degree  polynomials  as  will  be  seen. 
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For  the  viscosity,  we  assume 


so  that 

, JA  Z -0.36 

V IV  “ 

Combining  eq.  (22,  (23)  and  (25)  into  (20)  we  get: 


where : 


F, 

z foo 11  oo 

1 2lc‘-  ) C- 

' ( TruJ  Fl 

= 

f"(f. 

) ^oa  i 

(26) 

F = ( 

_ 

V-H 
' Z ) 

(- 

'A® 

P"  CD  / 

(27) 

i 

' Z I 

V 

- 0-38 

. 1 
l = 

b 4/y-.) 

y-» 

8 

* — ) 
II 

(28) 

The  value  of  F was  computed  for  = 2, 3, 4, 5,  and  6,  and  for  y = 1.2,  1.3, 
1.4,  1.5,  1.6  and  1.667  at  each  M,  and  is  shown  on  Table  III.  The  integral 
i was  computed  numerically  using  the  TSS  NUMINT  program,  with  M expressed  as 
a function  of  (^/^)  through  eq.  (12)  and  eq.  (17).  It  is  shown  on  Table  IV. 

It  is  important  to  relate  the  viscous  force  on  the  nozzle  with  identifiable 
properties  of  the  nozzle  flow.  From  the  Karman  integral  equation  for  the 
momentum  thickness  on  a plate,  we  get  one  approximate  relation: 


„ , 1 Jb 

r®u»  jj 


from  which, 


f * 

Zjo  Tu A]  -2?aul  Qx 
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where  9*  is  the  momentum  thickness  on  the  nozzle  (assumed  equivalent  to  a 
flat  plate)  trailing  edge.  Inserting  (30)  into  (26)  we  get 


A\  =7^-.  = ( -±b)n 

1 - poouoo  V ' 


(31) 


Another  relation,  to  prove  very  useful  later  on,  connects  Fv  with  the 
"Momentum  defect"  of  the  wake: 

(Viscous  force)/span  = (Momentum  wake  change)/span 

A 


or 


CD 


(<r.  tine  w A We) 


(32) 


too 


When  divided  by  l/2p  u^  , the  r.h.s.  is  the  familiar  momentum  defect  of  the 
wake : 


.on 


piA 


- (1-rWl 


(33) 


)0  ?oou^ 

from  which,  with  the  aid  of  (30),  (31)  and  (32): 


cjx  = 4Qt 


(34) 


Thus  we  find  that  the  (normalized)  viscous  force  of  the  nozzle  cusp  is  equal 
to  the  momentum  defect  of  its  wake  CDh  and  also  approximately  equal  to  four 
times  the  boundary  layer  momentum  thickness  at  the  nozzle  trailing  edge. 
Although  this  section  deals  with  laminar  flow  on  the  nozzle  cusp,  equations 
(30)  and  (34)  are  valid  for  turbulent  flow  as  well. 

As  a result  of  the  above,  we  can  now  combine  (31)  and  (34)  into 


-f  _ / 2.14.  \ 


feol 


)Fi 


(35) 


In  the  GDL  problem,  in  which  the  nozzle  axes  of  symmetry  are  separated  by  a 
distance  hi  according  to  Fig.  1,  it  will  be  convenient  to  refer  all  distances 
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to  hi  ; thus. 


CpURlflat 

21 C 


_ Cd^  l fot  _ j Cfcft  \ \/^®hi  / fi,  \ * p • 

' T;r  2cj  * > zcL  lr/ 


from  which: 


C-^ 'lR<ooh,  = 2C,.'/|Fi-=  JfMo,,!')  <37> 

with  2Cl®1.328.  This  last  equality  is  possible  because,  due  to  geometric 
similarity,  GDL  nozzles  have  a length-to-exit-height  ratio  dependent  on  M 
(see  Ref.  1). 

The  net  result,  expressed  by  (37),  then  is  that  the  momentum  defect  of  wakes 
shed  from  GDL  nozzle  cusps,  expressed  as  a fraction  of  the  nozzle  exit 
height  hi , is  inversely  proportional  to  the  square  root  of  the  exit  Reynolds 
number 

f?0  . = fooUoc'k.  (38) 

'-onto,  — “ “ 

V-oo 

and  directly  proportional  to  a function  of  the  exit  Mach  number  and  y. 
The  function  J has  been  computed  and  is  shown  on  Table  V.  The  inverse 
proportionality  to  Reynolds  number  is  fully  expected  for  laminar  flows.  As 
Table  V also  shows,  ? greatly  increases  with  Mach  number  and  generally 
decreases  with  y.  The  function  f is  also  shown  on  Fig.  3. 

6.2  HEAT  EXCHANGE  IN  LAMINAR  FLOW 

For  laminar  wakes  we  take  the  following  expression  to  represent  the  energy 
flux  (energy/area/time)  from  the  nozzle  surface: 

<39) 

The  total  energy  exchange  per  unit  span  will  then  be 

l , f‘ 

<*pfpT0j(i)  <4°> 

o ° 
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where 


(41) 


so  that  if  TW<TQ  (cooled  plate),  Q < o also.  Thus  Q is  negative  for  heat 
flowing  into  the  nozzle  and  positive  if  heat  is  flowing  from  the  nozzle  to 
the  fluid. 


From  (40)  we  get, 


for  c 

P 


constant , 


(42) 


Physically  this  parameter  is  the  ratio  of  the  heat  exchanged  between  nozzle 
cusp  and  flow,  to  the  heat  content  of  the  flow  exiting  the  nozzle.  For 
later  purposes,  it  is  necessary  to  know  the  corresponding  ratio  for  the 
kinetic  "power"  exiting  the  nozzle: 


Q-  / Q.  \ I "loo 

i Poou®  Ct> fcouoo(f"^»  ^oo  J 


( a . ) 

Cp  \ 

1 — — -«•  1 

1 Q.  > 

1 ' (43) 

fcoU®^p  "^0 

'2  coh 

If  the  integral  in  (42)  is  called 


l 


I 

f 

I 


r+i 

i/m  f.  V,"] 


> O 


i Vmo,,  y) 


(see  Table  VI  for  i'  values)  then 


(44) 
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(45) 


a = 


Using  (35),  we  get 


*Jcu  / ft 


y-fi 

'-38  +47p) 

(t  ) 


'Ll 


£2/3{*^l  V^Tco 


1.38 


Poo 


’<1^-0  • ' 


-rjr  = 5*  (M^y  Y) 


l F 


or 


<*  GV3  (y-i)M* 

Q.  ■=  « S - ( ^ -«)  SfMa,,  Y) 


(46) 


(47) 


The  quantity  Q,  which  will  feature  prominently  in  the  prediction  of  wake 
temperature,  thus  depends  on  a new  entity,  the  wall-to-stagnation  tempera- 
ture ratio,  and  the  parameter  S.  The  latter  can  be  considered  to  be  the 
"efficiency"  of  cooling  or  heating  the  wake  by  cooling  or  heating  the  nozzle; 
it  is  shown  on  Table  VII  and  Fig.  4.  As  with  Q,  Q<0  if  the  nozzle  is  cooled 
(TW<T0)  and  Q>0  if  the  nozzle  is  heated.  Note  that  Q = -S  is  the  minimum 
physically  attainable,  within  the  present  approximations. 


The  results  of  Fig.  4 and  Table  VII  indicate  that  S is  highest  at  low 
and  small  y,  and  lowest  at  high  and  large  y.  This  result  has  also  been 
verified  for  the  simple  case  of  a flat  plate  at  constant  flow  conditions 
U®,  etc.  Using  the  Reynolds  analogy 


r - £i  (48) 

2 <5 

for  laminar  flow  and  the  relation  between  the  T.E.  momentum  thickness  0*^ 
and  coh  (q*j  = Cjjh/4),  it  was  shown  that 

5:  _ !_.  r L 4.  _l'|  («) 

* S'  L (y-o  z 1 

which  has  the  same  qualitative  behavior  with  Mco  and  y as  the  curves  of  Fig.  4. 
This  result  is  physically  connected  to  the  definition  of  Q (see  Eq.  (42)  ) 
which  is  the  proper  definition  for  calculating  wake  properties  (Section  7). 

6.3  VISCOUS  DRAG  IN  TURBULENT  FLOW 

As  for  the  laminar  case,  the  total  viscous  drag  of  the  nozzle  cusp  per  unit 
span  is  given  by  Eq.  (5): 
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where  will  now  be  taken  from  Ref.  2 to  be  typical  for  turbulent  flow: 


S’ 

7 I 


ct=  we-;  = cTf( 


’ M)  • 


t + 1 ' 1 <*■ 

The  dependence  of  c,  on  turbulent  heat  transfer  is  neglected  consistently 
with  assumption  D,  Section  3.  If  Cf  is  then  re-written  as 


(f^oo  l) 


'h  u V7 


we  get 


^ = 2* k k b* (53> 


% ‘ ' f * ' ^ P-oo  i P0 


This  expression  can  be  further  reduced  to 


C*  - H * ^ ^ v (54, 

p»  < teo.(  /7'  a I \ 2 I ' 


where 


(t)  m P-  m -Mi) 


+ 1 U' 

This  integral  has  been  calculated  and  is  shown  on  Table  VII.  The  factor  of 
this  integral,  furthermore , can  be  also  reduced  so  that  (54)  finally  reads 

,/  )3  Y Y+l 

h 2 o 101 +—■ — 

'^V,  WL"  (34, 

2CT  l 1 

The  factor  N(Mro,Y)  appears  on  Table  IX. 

As  was  done  with  the  laminar  case  (eq.  (37)  ),  we  can  express  the  la:t  result 
in  terms  of  hi  as  follows: 


12 


Co*.(R«a>0 


L - 1C 


/»!ff 

T (t"!  Nl  l = G O^ao  > Y ) (57) 


with  2Ct  = 2 x 0.0262  = .05424.  This  quantity  is  tabulated  in  Table  X and 
also  plotted  on  Fig.  5.  It  is  interesting  to  compare  Fig.  5 with  Fig.  3 and, 
considering  Eq.  (34),  to  see  how  the  boundary  layer  momentum  thickness  Q*^ 
at  the  T.E.  behaves  in  laminar  and  turbulent  flows.  It  is  seen  that  Cph/fij 
(and  9*^)  in  both  cases  increases  with  and  decreases  with  y.  Quantita- 
tively, however,  there  are  differences  in  magnitude  between  the  two  cases 
because  of  the  Reynolds  number  exponent  (compare  Eq . (34)  with  (56)  ).  For 
example,  if  Re^hi  = 10^*: 


Laminar  case:  = 0.00&2. 

Turbulent  case:  = 0.0125 


(58) 

(59) 


However,  in  practical  circumstances  Re^h^  is  usually  much  larger  if  the  flow 
is  turbulent;  thus  for  a laminar  flow  with  Reoohj  = 10^  the  c^h/h  is  larger 
(=0.0196)  than  it  is  for  a turbulent  flow  with  Reco^  = 107  (=0.0088). 


6.4  HEAT  EXCHANGE  IN  TURBULENT  FLOW 


In  contrast  to  the  laminar  case,  the  heat  transfer  coefficient  in  the  turbu- 
lent case  obeys  the  Reynolds  analogy  directly: 


c = Si 
H 2. 


(60) 


Thus,  with  the  Prandtl  number  missing  in  the  heat  flux  equation  (40),  the 
dimensional  heat  transfer  rate  from  the  entire  nozzle  is 


a 


= 

Too  J 


pu 


foo^ao^p 

with  q/  defined  again  as  in  (41).  Thus, 

a ^ 


0 


c, 


-(I) 


(61) 


pu 


Cr  d. 


f) 


2 f(»  U00( 

Using  the  turbulent  skin  friction  coefficient  from  Eq.  (51),  we  obtain 

■ in 


(62) 


T • m ^ 

^ = li-  ^ ^ (Mao,  y) 
N t 


(63) 


oc 


o£ 


*These  examples  are  computed  for  Mra  = 4,  y = 1.4 
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where : 


3 

7 


(64) 
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H-o  ' 

1 (m) 

'ft*/  \ 

■F® j 

f \ X J 

r+i 

Y-' 


y-1  x 

n« 


M, 

t 3 jr+i 


- ■ KfM  V’ ^ ( jl 
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N from  Table  IX 


(65) 

(66) 
(67) 


i"  ^rom  Eq.  (54)  and  Table  VIII. 

The  quantity  Q /a,  which  is  again  considered  as  a heating  "efficiency"  of  the 
nozzle,  appears  in  Table  XI  and  Fig,  6.  As  in  the  laminar  flow  case  (Fig. 4), 
Q/q.  decreases  with  increasing  M^  and  also  with  increasing  y.  Comparing 
Fig.  4 with  6,  we  also  observe  that  the  heat  exchange  is  more  efficient  in 
turbulent  flow  but  only  at  the  higher  MOT  and  larger  y.  In  practical  terms, 
however,  caution  in  interpreting  this  result  is  needed  because  of  the 
manner  of  defining  Q (see  Eq,  (62)). 


7.  THEORY  OF  WAKE  DECAY 


7.1  LAMINAR  WAKES 

7.1.1  KUBOTA'S  LINEARIZED  SIMILARITY  THEORY  FOR  COMPRESSIBLE  WAKES 

Kubota  (Ref.  3)  has  discovered  similarity  solutions  of  the  boundary  layer 
equations  as  used  to  describe  two-dimensional,  laminar,  compressible  wakes. 
For  completeness,  his  results  are  given  here  in  outline.  We  will  adopt 
Kubota's  notation:  subscript  0 (as  in  ( )q) refers  to  the  wake  conditions  at 

some  arbitrary  initial  station;  ( )e  refers  to  conditions  just  external  to 
and  ( )„  to  conditions  very  far  away  from  the  wake.  An  overbar,  ( ) refers 
to  physical  (dimensional)  coordinates.* 

Kubota’s  coordinates  are: 


Longitudinal  distance: 

1 f * Pe  P'e  ue  . - 
x = — — _ J Y 

(68) 

Lateral  distance: 

a - u>  j & j5  f % 

(69) 

*This  subscript  notation  applies  to  this  Section  only.  In  Section  7.1.2  we 
return  to  our  own  notation. 
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where  d is  a characteristic  (constant)  length  and 


Re 

III 

1 8C| 
a- 

(70) 

hot 

His  solutions  are:  a 

ue-u  A _ 4* 

V A / — «.  Q 

S'  (*dUe  e 

(71) 

Ue  u|fx 

2uf  J0  Vex  +(i-er)£'  (? 

» T-Te  B'  -ClS/4>t 


He 


( i 


, -ejn*  A'  - , 

-I7,e  ) 


(72) 


(73) 


(74) 


(75) 


The  quantities  w,  h and  9 are  the  fluid  velocity  and  temperature  and  the 
momentum  thickness,  respectively. 


7.1.2  ADAPTATION  OF  KUBOTA'S  THEORY  TO  THE  PRESENT  PROBLEM 

In  this  section,  we  will  re-write  Kubota's  formulas  (68-75)  adapting  them  to 
our  nomenclature  and  subjecting  them  to  simplifications  suitable  to  the  GDL 
nozzle  problem. 

In  the  present  problem  we  shall  assume: 

- "The  nozzle  exhaust  flow  outside  the  wake  is  everywhere  uniform." 

This  assumption  approximates  well  the  flow  exiting  ideal,  well-designed 
nozzles  with  thin  trailing  edges  in  parallel  flow.  Thus,  quantities  sub- 
scribed by  "0"  or  "e"  in  eqs.  (68-75)  become  equal  to  those  subscribed  by  """ 
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Still  within  his  nomenclature  but  under  our  present  conditions,  Kubota's 
equations  (67)  and  (68)  then  become 


X 


(76) 


where  y”  is  the  familiar  (transformed)  compressible  coordinate: 
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'T 

)oo 

(77) 

so  that,  also 

cU,  = - 

f 

^ 'A 

(78) 

J s 

fee 

a 

Kubota's  central  approximation  ("linearization")  consists  of  assuming,  as 
we  will  also  in  this  work,  that 

w«l  (79) 

Then  his  momentum  thickness,  from  Eq.  (74)  is 
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(80) 


where,  from  Ref.  4,  c^h  is  the  "momentum  thickness"  of  the  wake: 


c0l  = 4 


0O  - - 

O fan 


acb 


(81) 


We  thus  see  that  Kubota's  "momentum  thickness"  q is  actually  equal  to  the 
momentum  thickness  g*,  of  the  boundary  layer  on  the  T.E.  of  a flat  plate 
generating  the  wake  (see  Eq.  (34)  also). 


Another  modification  to  Kubota's  theory  has  to  do  with  x.  In  our  nomenclature, 
Ax  is  the  distance  downstream  of  the  wake  "virtual  origin",  located  at  xQ. 
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(82) 


Thus,  Kubota's  x should  be  changed  to 

Ax  = x-xQ 


Because  of  previous  assumptions,  pressure  gradients  around  the  wake  are 
neglected: 


d_  = d_ 
dx  dy 


0 


(83) 


We  are  now  ready  to  obtain  the  flow  properties  in  the  wake.  Overbars  on 
physical  (dimensional)  properties  are  excluded  from  here  on.  Also,  from 
here  on  we  return  to  our  own  nomenclature  (see  list  of  symbols  and  Fig.  J_) . 


7. 1.2.1  VELOCITY  DEFECT  FOR  LAMINAR  FLOW 

Using  Eq.  (82)  and  (83),  in  (71)  we  obtain  on  the  wake  axis  (y  = y = 0): 


Thus, 


Uo,-U(0> 

< 

(84) 

^*00 

* f fi  - 

u-  57  'll 

(85) 

cA  f*£‘  J_ 

4 v TC.  VT* 

(86) 

where  Re'  = unit  Reynolds  number  = D u^/^,  In  the  GDL  problem  it  will  be 
convenient  to  normalize  distances  wit?  h;  , the  nozzle  exit  height  (Fig.  1). 
Thus , 


W 


4 ft,  \/tc 


A*  j 


(87) 


and  using  Eq.  (37),  we  find  that 


(88) 


Thus  the  velocity  defect  on  the  wake  centerplane  depends  only  on  the  down- 
stream distance  and  on  the  exit  Mach  number  (and  y).  Surprisingly,  there  is 
no  Reynolds  number  dependence.  If  the  virtual  origin  xQ  is  known,  w can 
thus  be  found  immediately  from  the  f values  of  Table  V,  as  computed  in 
Section  6.1. 
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7. 1.2. 2  TEMPERATURE  DEFECT  IN  LAMINAR  FLOW 

At  y = y1  = 0,  Eq.  (72)  gives,  combined  with  (88): 

^ _ T(0)— Top  g ' 

Too  \l  Cyc 


(89) 


with 


B 


a 


^ "2.  ^cq^co  ^p"Tao  V 


&Rt' 

TC 


(90) 


Thus, 

t • {£  o f-)Mi  i 


where 


a 


i^e 


(-  D~fxj 


_ ^qo  ^ cm 


^-oo 


(91) 


(92) 


Invoking  Eq.  (37)  and  (42),  we  finally  get 


(,4qO W <«) 

"Tqj 

Thus  the  temperature  on  the  wake  centerplane  is  seen  to  depend  on  the  laminar 
Prandtl  number  q,  ratio  of  specific  heats  y,  exit  Mach  number  and  the  non- 
dimensional  heat  exchange  rate  Q,  defined  in  Eq.  (42)  and  tabulated  in 
Table  VII  (see  also  Fig.  4).  Note  that,  from  Eq.  (47), 


a <0  IF  / To  < 


(94) 


Thus,  t decreases  if  the  nozzle  surface  temperature  decteases  below  T0. 

7. 1.2. 3 PRESSURE  DEFECT  IN  LAMINAR  FLOW 

Consistent  with  Kubota's  assumptions,  the  pressure  everywhere  inside  the  wake 
is  assumed  to  be  constant  and  equal  to  the  GDL  nozzle  exit  pressure  p^. 

7. 1.2. 4 DENSITY  DEFECT  IN  LAMINAR  FLOW 


Because  of  the  constant-pressure  assumption,  the  density  defect  can  be  com- 
puted directly  from  the  temperature  defect  of  Eq.  (93): 


)®  -t-H 


(95) 
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Thus,  the  density  defect  does  not  show  the  (Ax/hj )” * decay  behavior  exhibited 
by  the  velocity  and  temperature  defects. 

7. 1.2. 5 TOTAL  TEMPERATURE  DEFECT  IN  LAMINAR  FLOW 

From  Kubota’s  Eq,  (73),  using  our  nomenclature,  we  obtain  at  y = y = 0: 


To(0)  — "focD 

\TX 

/ / A V 

(96) 

~^0QD 

' 

■ (b  - — ) 

Too  ' 

where 

, again 

» + 

CD  / 2- 

(97) 

Inserting  the 

expression  for  Q 

(Eq. 

, (42)  ) 

into  B’  (Eq.  (73)  ),  we 

get 

/ 

B 

r , 
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(98) 

Also, 

f rom  Eq , 

. (73) 

for  A' : 

A/ 

■ (v- 

Cb~^  vj  ££ 

(99) 
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Since 

Re  5 P»l 
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, and  using 

Eq. 

(88)  for  w. 

& = 

1 

^oo 

. p" 

V Ay  4 

~ lY-OiAs 

(100) 

or 

<0  = fyj.(  l+a)_n  doi) 

where  w is  again  given  by  (88). 
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7. 1.2. 6 LATERAL  LAMINAR  VELOCITY  VARIATION,  AND  THE  TRANSVERSE  SCALE 
From  Kubota's  Eq.  (71)  we  get,  in  our  nomenclature 


- Ucn-U  w '4*  [ «?e 

u.  = — - = € = — 

-u(o)  L -t  d2- 


(102) 


Inserting  (86)  into  (102)  we  get: 


~ f “mi  \ 

u = 4 


(103) 


Define  a transverse  scale  L and  a dimensionless  lateral  coordinate  T)  by 

t.  ~ 

( L L io  * 


(104) 


Thus , 


-0.-78  Si 


u = e 


= e 


(105) 


Note  that  since  cph  is  a function  of  the  Reynolds  number  according  to  Eq.  (37) 
and  Fig.  3,  the  lateral  velocity  variation  is  also  a function  of  the  Reynolds 
number, 

7. 1.2. 7 LATERAL  TEMPERATURE  VARIATION  IN  LAMINAR  FLOW 

From  Eq.  (72)  we  arrive  at  the  counterpart  of  Eq.  (105)  for  the  temperature: 


T s = e 4 

TfcA  - Too 


_o.7gs;£'rj 


= e 


(106) 


7. 1.2. 8 LATERAL  PRESSURE  VARIATION  IN  LAMINAR  FLOW 


The  static  pressure  p^  is  constant  across  and  along  the  wake. 


7. 1.2. 9 LATERAL  DENSITY  VARIATION  IN  LAMINAR  FLOW 


Because  of  the  constancy  of  pressure,  the  lateral  density  variation  is  given 
by  ^ 

f«-f  T(-t-H) 


T-fc  + 


(107) 


with  T and  t supplied  by  Eq.  (106)  and  (93),  respectively.  Thus,  as  is  well 
known  (Ref.  4)  the  lateral  density  profile  is  not  similar  (i.e.,  ■jy  depends 
on  /\x  through  t)  until  t = 0,  which  of  course  occurs  only  asymptotically, 
far  downstream. 

7.1.2.10  LATERAL  TOTAL  TEMPERATURE  VARIATION  IN  LAMINAR  FLOW 
From  the  formulas  presented  so  far,  we  can  show  that  since 


**  -t 

T _ o 'Ooo 
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(108) 


then 


~ !4-a)e  4 1 - e 4 N 

T0  = — =r  T0(ri  ; S',  0.) 

Vff  ( I + GL)  - I 


(109) 


which  is  a similarity  solution. 

Fig.  7 shows  the  lateral  profiles  of  velocity,  temperature,  density  and  total 
temperature  In  the  laminar-wake  case,  the  last  three  for  <j  = 0.75. 

7. 1.2. '11  APPROXIMATE  CALCULATION  OF  LAMINAR  WAKE  THICKNESS 

We  define  the  "transformed  velocity  half  thickness"  as  the  lateral  distance 
^u  at  which  IT,  from  Eq.  (105),  has  decreased  to  a certain  level.  The  choice 
of  this  level  is  arbitrary,  and  will  be  taken  here  as  0.01: 


IA  = O.OI 


Then  from  Eq.  (105)  it  can  be  shown  that 

A/ 

r 

6 


^ =4.z<n) 


where 


KV  - 


i 1 1 

00  ‘uo 


no 


(110) 


(111) 


(112) 
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Because  of  symmetry  about  the  wake  centerplane  (axis),  then  the  (full) 
velocity  thickness  is 

— E Z = 8.  S8  (113) 

t,  €,  N ReEh. 

with  analogous  results  for  the  transformed  temperature  thickness: 

hi  - tfl  \fhil l (H4) 

^ V *6^, 

It  will  be  desirable  to  have  a formula  for  the  physical  thicknesses.  This 
is  impossible  to  do  in  the  context  of  similarity  because  the  density  ratio 
depends  on  x (see  Eq.  (107)  ). 

A serviceable  approximation  can  be  had  by  assuming  that  this  ratio  varies 
linearly  with  y: 


iZ  - + T. 


where  yfc  is  the  physical  wake  "edge."  Then, 
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and  since,  from  (107) 


fJO)  _ t j-  2. 

f oo  ‘ 


we  get  for  the  plysical  thicknesses, 


U/i,  = \J~EZL. 


(119) 


— T ^ 

*,  = **  *, 


(120) 


These  quantities  are  the  only  ones  calculated  so  far  for  the  laminar  wake, 
which  depend  on  the  Reynolds  number. 


7.2  TURBULENT  WAKES 

7.2.1  LINEARIZED  SIMILARITY  THEORY  FOR  TURBULENT  WAKES 


Since  no  straightforward  turbulence  theory  is  possible,  there  exists  no 
formal  counterpart  of  Kubota's  theory  which  can  be  developed  from  first 
principles  when  the  wake  is  turbulent.  There  exists,  however,  a linearized 
similarity  solution  for  turbulent  flows  for  the  "far  wake"  limit  of  w «1 
(Ref.  4).  This  solution  involves  a single  assumption  which,  as  can  be  ex- 
pected, indirectly  concerns  the  turbulent  viscosity.  The  solution,  which  is 
supported  strongly  by  experiments  (Refs.  4,5,  and  6),  is  very  appropriate  to 
use  here  because  of  its  kinship  to  the  Kubota  theory  (through  the  w « 1 
approximation)  and  because  it  yields  formulas  almost  exactly  the  same  as 
presented  in  Section  7.1. 

The  starting  point  of  the  "theory"  is  the  same  as  for  Kubota's  theory,  e.g., 
steady,  uniform  isobaric  flow,  etc.  Then,  the  key  assumption  regarding  the 
velocity  defect  is  made: 


W 


where  Rx  is  the  so-called  turbulent  Reynolds  number. 
Assumption  (121)  now  leads  to  the  following  results: 
7. 2. 1.1  VELOCITY  DEFECT  IN  TURBULENT  FLOW 


(121) 


The  velocity  defect  as  given  by  (121)  can  be  further  expressed  in  the  desired 
coordinates  by  utilizing  Eq.  (57): 
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Thus , 


(123) 
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If  we  adopt  the  experimentally  known  value  for  R^  from  Ref.  4: 


Rt  = 13 


and  noting  that  C^,  = 0.0262,  we  obtain 

*,  \/Z 
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Vj 


where 
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(125) 


(126) 


The  quantity  K has  been  tabulated  in  Table  XIII  for  Re^  = 10^.  Note  that 
for  other  values  of  Re^  the  K will  not  change  much,  sihce  it  depends  on 
the  former  to  the  l/14th  power,  according  to  Eq.  (126). 


7. 2. 1.2  TEMPERATURE  DEFECT  IN  TURBULENT  FLOW 


t=  =V5t  (v->)mIo+gDw 

loo 

where  Q is  defined  exactly  as  before  (Eq.  (42)  ). 

7. 2. 1.3  DENSITY  AND  PRESSURE  DEFECTS  IN  TURBULENT  FLOW 

As  for  the  laminar  flow,  the  pressure  is  everywhere  constant,  and 

x_  P«rf“a  _ /t_. 

foo  t + l 

7. 2. 1.4  TOTAL  TEMPERATURE  DEFECT  IN  TURBULENT  FLOW 

z. 
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7. 2. 1.5  LATERAL  TURBULENT  VELOCITY  VARIATION 
From  Ref.  4 we  get 

~ -O.M'n2’ 
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(127) 
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(129) 
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7. 2. 1.6  LATERAL  TURBULENT  TEMPERATURE  VARIATION 


From  Ref.  4 we  also  get 


t = e 


-0.  4 -S-Y| 


(131) 


7. 2. 1.7  LATERAL  DENSITY  AND  PRESSURE  VARIATIONS 

Again,  the  pressure  is  everywhere  constant.  The  lateral  density  variation 
o'  is  derived  exactly  as  for  the  laminar  wake: 


(132) 


r f«-f(<»  Ti  + I 

7. 2. 1.8  LATERAL  TURBULENT  TOTAL  TEMPERATURE  VARIATION 
The  lateral  variation  of  the  total  temperature  is  again  given  by 

T — 1 -*-Qu)  T - (X.  __  To  -T0Cp  (133) 

( I 4-0.)  — ( T^(0)-Toao 

, ~ and  u are  the  turbulent  lateral  variations  (Eqs.  (130)  and  (131)), 

where  the  i 

A note  of  caution  £s  needed  in  regard  to  Eq.  (133)  (the  same  applies  to  the 
laminar  later  variation,  Eq.  (109)  ). 


It  can  be  seen  that  for 


Gl  — — 

fer 


(134) 


The  denominator  is  zero;  since  is  usually  of  order  0.7  (in  laminar  or 
turbulent  flow),  this  condition  occurs  at  Q ~ 0.2,  i.e.,  when  the  wake  (or 
nozzle)  is  slightly  heated.  Physically,  the  explanation  is  simple:  when 

heating  is  applied  the  wake  axis,  which  is  normally  at  T(0)  < T0  (as  can 
be  verified  from  (129)  when  Q = 0),  begins  rising  in  temperature  until 
To(0)  = To,  or  Tq(0)  - T.  = 0. 

7. 2. 1.9  APPROXIMATE  CALCULATION  OF  THE  TURBULENT  WAKE  THICKNESS 

From  the  experimental  results  of  Ref.  4,  we  note  that  the  transformed  "half- 
thickness" of  the  turbulent  wake  (defined  as  the  point  where  ij  = 0,01)  is 


^ K V\ 
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(135) 
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Utilizing  the  approximate  inverse  transformation  (117)  between  and  y , 
and  again  calling  the  "full  thickness"  bu  = 2yu,  we  get 


2.58 


■fc  + 1 

V+2 


(136) 


Then,  utilizing  the  relation  (125)  for  w in  turbulent  nozzle  flow,  we  get 
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The  quantity  K is  tabulated  in  Table  XIII  for  Reccj^  = 10^.  Note  that  for 
other  values  of  Rerah  the  K will  not  change  much  since  it  depends  on  the 
former  to  the  l/14th  power,  according  to  Eq.  (126). 

7.3  THE  VIRTUAL  ORIGIN  ASSUMPTION 


The  formulas  presented  so  far  express  the  variables  in  terms  of  the  longitud- 
inal coordinate  Ax  = x-xQ,  where  xQ  is  the  "virtual  origin"  of  the  wake. 

The  location  of  the  virtual  origin  cannot  be  found  from  first  principles, 
and  must  be  obtained  from  experiment.  In  the  meantime,  an  "educated  guess" 
of  xQ  is  in  order  to  that  some  numerical  applications  of  the  previous  results 
can  be  obtained. 


Since  the  wake  thickness  is  finite  at  the  cusp  T.E.,  it  is  certainly  incorrect 
to  locate  x0  there,  i.e.,  to  put  xQ  = 0 (see  Fig.  1 for  the  coordinate  con- 
vention). Since  the  nozzle  has  been  assumed  to  be  a flat  plate  for  con- 
venience, and  since  its  wetted  length  (for  friction  and  heat  transfer  cal- 
culation purposes,  see  Section  6)  has  been  assumed  equal  to  -t,  we  shall 
assume  that 

= - i (138) 

so  that  the  virtual  origin  for  both  the  laminar  and  the  turbulent  wake  is 
assumed  to  lie  at  the  location  of  the  nozzle  throat. 


The  above  assumption  means  that  as  it  appears  in  the  various  formulas,  the 
coordinate 
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so  that,  for  example,  a point  lying  one  unit  of  h;  (x  = h:  ) downstream  of 
the  T.E.  actually  lies  1 + downstream  of  the  origin.  Since  t/hj  is  also 
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a function  of  M“,  the  trailing  edge  itself  (x  = 0,  & x/hi = i/h: ) lies  downstream 
of  the  origin  by  an  amount  dependent  on  given  on  Table  XIV. 

7.4  THE  TRAILING  EDGE  NON- EQUILIBRIUM  WAKE  REGION 

It  is  clear  from  the  above  discussion  of  the  virtual  origin  that,  according 
to  the  formulas  given,  the  wake  will  have  already  diffused  at  the  T.E. 
position;  for  example,  for  laminar  flow  M^,  = 2,  y = 1*6,  w at  the  trailing 
edge  is  0.37.  Simple  physical  arguments,  however,  dictate  that  at  the 
trailing  edge  w = 1.  There  will  be  a region,  therefore,  in  which  the  wake 
adjusts  from  w * 1 at  the  T.E.  to  the  values  given  by  the  present  theory, 
and  beyond  which  w follows  the  Ax/hj  dependence  given  herein.  This  region 
is  called  the  "trailing  edge  (T.E.)  non-equilibrium  region."* 

The  extent  of  the  non-equilibrium  region  can  be  calculated  with  difficulty 
for  laminar  flows;  in  turbulent  flow  the  usual  hardships  about  turbulent 
viscosity  are  added.  In  either  case,  its  extent  is  of  course,  dependent 
on  the  size  and  shape  of  the  cusp  T.E.  No  such  calculations  have  been 
attempted  herein;  however,  since  the  extent  of  the  region  should  scale 
with  the  trailing  edge  height,  the  assumption  of  "sharp"  trailing  edge 
adopted  here  should  help  keep  the  non-equilibrium  region  short. 

It  should  also  be  stressed  that  problems  arising  because  of  the  non-equilibrium 
region  are  quite  distinct  from  the  problem  of  meeting  the  "linearization" 
criterion  w « 1.  That  is,  even  after  the  flow  "loses  memory"  of  the  T.E., 
the  w may  not  be  low  enough  to  satisfy  the  validity  criteria  of  the  theory 
by  which  it  was  computed. 


*In  the  extensive  wake  literature,  it  is  also  called  the  "near  wake  region 
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8.  APPLICATIONS  OF  THE  THEORIES  FOR  LAMINAR  OR  TURBULENT  WAKES 
8.1  PROCEDURES  FOR  COMPUTIN G THE  WAKE  CHARACTERISTICS 

To  compute  wake  decay  for  GDL  nozzle  cusps,  a set  of  inputs  is  necessary: 


o Ratio  of  specific  heat  y (140) 
o Exit  Mach  number  Moo  (141) 
o Exit  Reynolds  number  Re^  (=pcoUa>/u<x>)  (142) 
o Stagnation  (supply)  temperature  To  (143) 
o Nozzle  wall  temperature  (average)  Tw  (144) 
o Nozzle  exit  height  tv  (145) 


Note  that  Re«,  can  be  calculated  for  any  gas  if  To,  Moo  and  the  supply  pressure 
Po  are  known  instead.  Note,  also,  that  y,  M^  and  hi  alone  are  sufficient  to 
compute  the  velocity  defect  w,  since  the  latter  does  not  depend*  on  Rei,  Tw  or 
To. 

The  first  thing  to  do  after  the  parameters  (139)-(144)  are  known  is  to  de- 
termine whether  the  wake  is  laminar  or  turbulent,  since  the  results  will  be 
markedly  different  for  these  two  types  of  wakes.  Wake  transition  correlations 
are  available  for  2-dimensional  supersonic  wakes,  in  Refs.  7,  8 and  elsewhere. 
The  present  formulas  are  designed  to  handle  "wholly  laminar"  wakes  (laminar 
boundary  layer  on  the  nozzle  followed  by  laminar  wake)  or  "wholly  tu.'bulent" 
wakes  (turbulent  nozzle  boundary  layer  followed  by  turb  : ent  wakes.)  As 
already  hinted  earlier,  the  transitional  case  is  complex  and  will  be  addressed 
in  a later  report. 

Once  the  state  of  the  wake  is  thus  determined,  the  equations  of  Section  7.1 
(together  with  the  appropriate  tables)  must  be  used  for  the  laminar  case, 
and  those  of  Section  7.2  for  the  turbulent  case.  If  laminar  wakes  a/e  in- 
volved, it  is  necessary  to  first  compute  the  function  ^ (Table  V)  for  the 
given  values  of  M<»  and  V.  This  Table  already  tabulates  Sr  for  M^  = 2, 3, 4, 5 
and  6 and  for  Y = 1.2,  1.3  ....  1.6,  1.667.  (Linear  interpolation  is 
adequate  for  intermediate  values  of  M^  and  y.)  With  J-  thus  known,  w can 
be  computed  as  a function  of  (^x/hi ) from  Eq . (88);  note  the  definition  of 
/\x  and  hj  from  Eq.  (139)  and  Fig.  1.  The  value  of  w so  found  can  now  be 
used  into  Eq . (93)  to  give  t(^x/h]).  In  this  case,  the  laminar  flow  Prandtl 

*In  the  laminar  case,  the  Re^  dependence  of  Cph  cancels  the  ReA  dependence 
of  w.  In  the  turbulent  case  w depends  so  weakly  on  Re<t  (see  Eq.  (124)  ) that 
computation  at  a typical  Rei  will  suffice,  within  a small  percentage  error, 
fur  a wide  range  of  Re<*>. 
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number  u is  needed  (typically  ~ 0.7)  and  Q can  be  found  from  the  Q/,-y  of 
Table  VII  as  follows: 


Gt  = 


& ( ^ 
<*  v T0 


(146) 


where  Tw  and  To  are  assumed  given,  as  stated  above.  With  Mro,  y,  w,  L and  Q 
available,  the  density  and  total  temperature  defects  r and  C)  (Eqs.  (95)  and 
(96)  can  also  be  computed.  Note  that  the  pressure  within  and  around  the 
wake  is  constant  (and  equal  to  the  nozzle  exit  pressure)  everywhere  in  this 
formulation. 


Since  t has  been  computed  from  (93),  the  wake  thickness  can  now  be  found 
from  (119).  In  this  computation,  the  nozzle  exit  Reynolds  number  Re^  is 
also  utilized. 


The  computation  of  the  lateral  distributions  is  somewhat  more  involved. 
Having  computed  w,  we  next  need  Cj)h  for  computing  L from  Eq.  (104).  By 
Eq.  (37), 


V7_ 

V 


(138) 


With  L thus  known  from  (104),  we  can  compute  Tl,  from  Eq.  (104),  for  any 
(transformed)  distance  y off  the  wake  axis;  then  the  velocity  at  ^ is  com- 
puted from  (105).  The  same  procedure  is  followed  for  the  lateral  distribu- 
tions of  temperature  (from  (106)),  density  (from  (107))  and  total  temperature 
(from  (109)). 


If  the  wake  is  turbulent  the  new  pertinent  function  (instead  of  9-  ),  is  K, 
given  in  Table  XIII.  Then  the  required  properties  are  obtained  from  the 
formulas  of  Section  7.2. 


8.2  A NUMERICAL  EXAMPLE 


To  illustrate  the  use  of  the  formulas,  we  shall  compute  the  following 
example.  A set  of  GDL  nozzles,  of  exit  height  1 cm  and  kept  at  temperature 
500°K,  discharges  helium  gas  at  exit  Mach  number  3.  If  the  supply  chamber 
upstream  of  the  nozzles  is  kept  at  a pressure  of  100  psia  and  temperature 
of  1000°K,  let  us  compute  the  gas  velocity,  pressure,  temperature,  etc.,  at 
a point  on  the  wake  axis  of  symmetry  lying  at  4 cm.  downstream  of  the  nozzle 
cusp  trailing  edge  (for  laminar  as  well  as  turbulent  flow.) 

Solution 


The  quantities  given  are: 

(147) 

(148) 


Meo  = 3 

y = 1.667 
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P0  = 100  psia  (14? 

T0  = 1000°K  (15( 

Tw  = 500°K  (151 

hj  5 1 cm  (15* 

x = 4 cm  (151 

We  can  first  compute  the  nozzle  exit  Reynolds  number  Re^  from  standard 
references,  such  as  Ref.  9,  from  the  given  quantities  (147),  (148),  (149) 
and  (150): 


Reh  = 127,000 


(154) 


Also,  it  is  convenient  to  compute  at  the  outset,  the  quantity  /\x/hj  from 
Eq.  (139)  the  quantity  (147)  and  Table  XIV: 

^ = * t 2.05- 4.  <155) 

We  can  now  solve  the  problem  first  by  assuming  a laminar  wake  (which  is  in 
anyway  appropriate  since  Re^  is  small).  From  Eq.  (88),  and  using  Table  V, 
for  the  inputs  (147)  and  (148),  we  get  for  the  velocity  defect. 


W = _!_  (fe.QS')  Z(  4.  4S4  ) =■  O ,2-SS 
4 Vr r 


(156) 


To  compute  the  temperature  defect  from  (93)  we  first  compute  Q as  follows: 
inputing  (150)  and  (151)  into  Eq.  (41): 

a =-0.5  (157) 

while  for  (147)  and  (148)  Table  VII  gives 


= 0.4-feM  , Q.  - —0-2-3 


(158) 


Inserting  CT  = 0.75  (the  typical  laminar  Prandtl  number)  as  well  as  (148), 
(147), (156)  and  (158)  into  Eq.  (93)  we  get 

t = 1.02  (159) 

The  density  defect  then  follows  directly  from  (95): 

r = 0.505  (160) 

The  total  temperature  defect  9 can  also  be  found  from  (101) 


q = -0.128 


(161) 
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To  compute  the  wake  thickness  we  utilize  (159),  (155),  (152)  and  (154)  into 
Eq.  (119): 

bu  = 0.0793  cm  (162) 


The  physical  (dimensional)  quantities  of  interest  can  now  also  be  computed 
directly  from  the  last  few  formulas  given  above,  utilizing  the  definitions 
(84),  (93),  etc.  For  example,  for  the  velocity  u(0)  and  temperature  T(0) 
on  the  wake  axis : 

u(0)  = u00(l-w)  (163) 

T(0)  = T^d  +t)  (164) 


The  stream  quantities  u„  and  T,*,  are  first  evaluated  for  isentropic  flow 
(using  the  given  quantities)  and  Ref.  9. 


Uoo  = 2.79x10^  cm/sec 

(165) 

T»  = 250°K 

(166) 

Then, 

u(0)  = 2.08x10'’  cm/sec 

(167) 

T(0)  = 505°K 

(168) 

For  the  turbulent  wake,  the  velocity  defect  is  obtained  from  inserting  the 
value  of  K from  Table  XIII,  and  input  (155),  into  Eq.  (125): 

w = 0.0563  (169) 


where  the  present  value  of  Rejlj  (input  (154))  and  that  used  in  setting  up 
Table  XIII  (Rehi  = 10^)  are  close  enough  to  make  practically  no  difference. 


For  the  temperature  defect,  use  is  made  of  the  same  a as  before  (see  (157) 
and  Q/y  from  Table  XI: 

Q = 0.5  (0.714)  = -0.357  (170) 


Thus  t is,  using  CTx  = 0.7,  from  Eq.  (127): 

t = \|o3  ^ (0.&4-3U 0.0S43)  =0.  I0Z 

Then  r from  Eq.  (127)  is 


O.l? 2 

1+0.132. 


= O.  IS*+ 


while  f)  from  Eq.  (129)  is 

P)  = -0.039 

i 
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(171) 


(172) 


(173) 


■ 

! 


J 


For  the  turbulent  wake  thickness  we  utilize  (171),  (152),  (155)  and  the 
value  of  K from  Table  (XIII)  into  Eq.  (137): 


b =.  0,Sl~7  - = 0.364  Cm.  (174) 

u 2 > I 6 2. 

Dimensional  examples  of  variables  are: 

u(0)  = 2.64xl05  cm/sec  (175) 

T(0)  « 295°K  (176) 

The  results  for  the  laminar  and  turbulent  wake  can  then  be  summarized  as 
follows : 


Quantity 

Laminar 

Turbulent 

Uoo-u(O) 

Velocity  defect  w = 

uoo 

0.255 

0.0563 

Velocity  u(0)  (cm/sec) 

2.08xl05 

2.64xl05 

Temperature  defect  t = (T(0)-Tao)/T0o 

1.02 

0.182 

Temperature  T(0) 

505 

296 

Density  defect  r = (pra-p(0)  )/pa3 

0.505 

0.154 

Total  temperature  defect  fi)  = (To(0)-To  VToa, 

-0.128 

-0.039 

Thickness  (cm) 

0.0793 

0.366 

By  their  definition  the  "defects"  shown  (of  velocity,  temperature,  etc.) 
are  measures  of  non-uniformities  in  the  GDL  cavity;  physically  they  resemble 
"troughs"  cutting  through  the  cavity.  We  see  that  the  laminar  case  is  dis- 
tinguished by  large  defects  and  thin  wakes;  the  inverse  is  true  for  the 
turbulent  case.  Thus  when  the  wakes  are  laminar  the  troughs  are  narrow  and 
deep,  while  in  the  turbulent  case  they  are  broad  and  shallow. 

8.3  NUMERICAL  CALCULATIONS  FOR  Moo-4,  y»  1.4 

Detailed  numerical  computations  for  the  laminar  and  turbulent  cases  have 
been  done  for  the  case  Moo  = 4,  y = 1«4  and  Re^  = 10^,  using  standard  values 
for  the  Prandtl  number  (cr  = 0.75,  = 0*7).  These  parameters  were  chosen 

to  correspond  to  those  of  current  experiments  in  the  Aeronutronic  Ford  Super- 
sonic Wind  Tunnel. 

The  results  of  this  calculation  are  shown  on  Figs.  8 through  12.  Fig.  8 
shows  the  velocity  defect  w for  the  laminar  and  turbulent  case;  this  quantity 
is  independent  of  Tw  for  both  cases.  There  is  no  Reynolds  number  dependence 
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in  the  laminar  case,  while  in  the  turbulent  case  the  Re^ -dependence  is  ex- 
tremely weak.  Note  the  large  difference  (by  a factor  of  5)  between  the  two 
cases.  In  the  turbulent  case,  too,  at  a downstream  distance  of  10  nozzle 
exit  heights  from  the  virtual  origin  (7.5  from  the  T.E.)  the  velocity  on  the 
wake  axis  differs  from  the  stream  velocity  by  only  a few  percent. 

The  temperature  defect  t is  shown  on  Fig.  9.  Like  w,  t decays  like  the 
inverse  square  root  of  distance  (obscured  in  the  Figure  because  of  the  semi- 
log plot).  Like  w,  t is  only  weakly  dependent  on  Re^  (and  then  only  in  the 
turbulent  case.)  Note  that,  as  expected,  cooling  the  nozzles  also  cools  the 
wake  in  almost  an  exact  proportion  to  the  nozzle  temperature  drop.  In 
absolute  magnitude,  however,  the  wake  cools  down  less:  for  example,  when 

the  nozzle  temperature  is  lowered  360°R  below  a supply  temperature  of  600°R, 
the  turbulent  wake  axis  temperature  at  Ax/hi  = 10  is  lowered  by  only  30°R. 

Fig.  10  shows  the  density  defect.  As  can  be  seen  from  Eqs.  (95)  and  (128) 
the  density  does  not  obey  the  inverse-square-root  decay  law  with  (Ax/hj). 
Also,  because  of  the  functional  form  of  the  Eqs.  (127)  and  (128),  the  de- 
pendence of  r on  Rehi  is  now  practically  nil  in  the  turbulent  case. 

It  should  be  stressed,  in  regard  to  Figs.  8 through  10,  that  the  various 
"defects"  plotted  are,  as  defined,  measures  of  the  flow  non-uniformity  in 
the  cavity  (and,  as  such,  useful  concepts  for  GDL  problems).  Thus,  their 
general  decrease  as  Tw  decreases  below  To  does  not  mean  that  the  temperature 
(or  density)  decreases,  but  rather  that  the  non-uniformities  in  the  GDL 
cavity  decrease  as  the  nozzle  temperature  decreases  in  relation  to  the 
supply  temperature.  Physically  this  happens  because,  while  the  main  gas 
flow  is  cooled  by  the  expansion  process,  the  wake  fluid  retains  the  friction- 
al heat  acquired  on  the  nozzle  surfaces  and  tends  to  be  much  hotter  than  the 
surrounding  fluid.  This  can  then  be  offset  by  keeping  the  nozzles  cooler 
relative  to  the  supply  chamber.* 

One  type  of  defect  which  increases  as  Tw  is  lowered  is  that  of  the  total 
temperature,  plotted  on  Fig.  11.  We  know  from  Eq.  (129)  thau  fi)  would  be 
identically  zero  for  unity  Prandtl  number  when  Tw  = To.  The  departure  of 
q-  ( and  (jj)  from  unity  creates  a finite  fi),  as  shown,  even  for  Tw  = To.  As 
the  temperature  Tw  is  lowered  p)  increases  negatively,  consistent  with 
energy  considerations. 

Wake  thicknesses,  as  defined  in  the  previous  formulas  (see  Eq.  (137)  for 
example,  are  plotted  on  Fig,  12.  Note  that  the  turbulent  wake,  irrespective 
*From  the  definition  of  t and  r (e.g.  Eq.  (127)  and  (128)  it  is  clear  that  t 
can  attain,  in  the  limit,  a minimum  value  of  -i.  When  near  this  limit,  t 
can  cause  r to  be  very  large  and  negative  in  sign,  implying  p(0)>>p0O.  This 
is  physically  realistic  and  implies  a very  dense  wake  in  this  case. 
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of  Reynolds  number,  is  about  6 times  thicker  than  the  laminar  wake  computed 
for  Re^  = 10^,  which  is  quite  possibly  at  the  lowest  end  of  the  GDL  range- 
of-interest.  The  laminar  wake  depends  on  the  inverse  square  root  of  Re^  , 
so  that  any  increase  in  Rej^  will  tend  to  "thin  down"  the  laminar  wake.  In 
either  laminar  or  turbulent  case  it  should  be  kept  in  mind  that  by/hj  does 
not  grow  as  the  square  root  of  (Ax/hj),  because  of  the  inverse  transformation 
(see  Section  7.1.2.11). 

By  a distance  equal  to  10(Ax/h1  ) it  is  seen  by  Fig.  12  that  the  laminar  wake 
occupies  about  10%  of  the  width  of  the  GDL  cavity  (by/hj  — 0.1),  but  the 
turbulent  wake  occupies  about  607.  of  the  cavity.  It  is  clear,  therefore, 
that  the  turbulent  wake  could,  in  some  instances,  totally  "fill"  the  GDL 
cavity. 

9 . CONCLUSIONS 

The  device  of  a "universal  family"  of  DeLaval  nozzles  has  enabled  us  to 
connect  the  integral  nozzle  properties  with  the  wake  properties.  Linearized 
similarity  formulas  were  then  derived  by  which  the  important  wake  quantities 
can  be  found  if  the  GDL  "design"  parameters  are  known.  Linearization  implies, 
in  practical  terms,  that  the  solutions  are  valid  only  when  w is  small,  that 
is  not  too  near  the  trailing  edge.  This  is  so  even  if  the  T.E.  is  ideally 
"sharp."  If  the  T.E.  is  truncated,  an  additional  reason  will  exist  for 
invalidating  these  solutions  immediately  downstream  of  it. 

The  most  remarkable  result  obtained  is  that  the  solutions  for  most  of  the 
laminar  wake  properties  are  Reynolds-number-independent.  This  "Reynolds 
number  independence  principle"  seems  to  have  been  first  discovered  by  Dewey 
(Ref.  10)  for  laminar  shear  layers  in  separated  flow  regions.  In  the  present 
case  it  arises  because  the  Re-dependence  of  the  nozzle  flow  cancels  the  Re- 
dependence of  the  wake  flow.  For  the  turbulent  case  there  is  a corresponding 
Reynolds-number  independence,  in  the  sense  that  once  a computation  is  done 
for  a typical  Re^  the  results  are  valid  for  other  Re^  within  a broad  range, 
to  good  accuracy. 

If  the  wake  defects  are  viewed  as  undesirable  inhomogeneities,  then  one  is 
better  off  to  design  for  low  and  large  y,  and  to  lower  Tw  considerably 

below  To.  The  Reynolds  number,  according  to  the  present  results,  is  of 
concern  only  when  the  wake  is  laminar,  and  even  then  it  affects  only  the 
wake  thickness. 

Finally,  it  must  be  stressed  that  when  the  wake  is  turbulent,  there  are 
fluctuations  superposed  on  the  mean  flow  field  as  calculated  herein.  These 
fluctuations  may  be  very  important  from  the  standpoint  of  optics,  and  will 
be  treated  in  a future  report. 
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TABLE  I 


TABLE  OF  PHYSICAL  COORDINATES  OF  CUSP  NOZZLES 


5 (in*> 

X(in*) 

? (in.) 

x(in.: 

O 

.5142 

.0398 

.4896 

.1808 

.4089 

.052 

.482 

.1916 

.4034 

.0575 

.4784 

.2029 

.3978 

.0618 

.4759 

.2153 

.3918 

.0646 

.4741 

.2272 

.3862 

.0674 

.4726 

.2398 

.3803 

.0699 

.4711 

.2530 

.3744 

.0734 

.469 

.2658 

.3686 

.0836 

.4627 

.2945 

.3563 

.0873 

.4603 

.3251 

.3437 

.0912 

.4579 

.3582 

.3307 

.0943 

.4562 

.394 

.3172 

.0971 

.4545 

.4326 

.3034 

.100 

.4527 

.4748 

.2889 

.1025 

.4511 

.5215 

.2742 

.1074 

.449 

.5715 

.2592 

.1125 

.4463 

.6274 

.2433 

.117 

.4436 

.6866 

.2274 

.1214 

.4407 

.7525 

.2111 

.1268 

.4376 

.904 

.1768 

.137 

.432 

1.0888 

.141 

.146 

.427 

1.3161 

.1045 

.1529 

.4235 

1.5902 

.0704 

.1616 

.4189 

1.9317 

.0407 

.1704 

.4143 

2.3479 

.0176 

2.6157 

.0076 

2.8545 

.005 
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TABLE  II 
COEFFICIENTS  CR 


m 

4 

i 

5 

6 

7 

1 

8 

i 

7.854531 

10.62409 

13.90497 

16.69116 

19.3553 

2 

-24.88887 

-51.14484 

-98.28814 

-156.8195 

-233.4886 

3 

31.898 

112.3132 

335.6014 

750.3165 

1499.482 

4 

-13.97012 

-111.6295 

-572.6051 

-1910.121 

-5325.4 

5 

_ 

40.77446 

471.7772 

2632.214 

10866.19 

6 

- 

- 

-149.4286 

-1851.714 

-12672.24 

7 

- 

- 

- 

520.4112 

7838.28 

8 

- 

-1991.188 

TABLE  III 


Values  of  F (Mco,  y)  for  Various  Values  of  M,,,,  y 
(Laminar  Viscous  Force  Calculation) 


SB 

2 

3 

4 

HSR 

6 

1.2 

1.199222 

1.977851 

4.011097 

17.32144 

1.3 

1.316166 

2.099284 

3.830711 

6.982828 

12.22896 

1.4 

1.4322 

2.231281 

3.804024 

6.361286 

10.17771 

1.5 

1.547565 

2.369352 

3.852264 

6.082258 

9.178528 

1.6 

1.662436 

2.511236 

3.942517 

5.974535 

8.654355 

1.667 

— 

1.739188 

2.607782 

4.013111 

5.959569 

8.451732 
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TABLE  IV 


Values  of  the  Integral  i for  Various  Values  of 
Mco  and  y (Laminar  Shear  Force  on  Nozzle) 


TABLE  V 

Values  of  the  Non-Dimensional  Parameter^  (M^,  y) 
(Laminar  Viscous  Force  Calculation) 


2 

3 

3.124 

4.653 

3.148 

4.579 

3.169 

4.528 

3.188 

4.495 

3.205 

4.47 

3 -214 

4.456 

7.092 


6.532 


6.202 


5.987 


5.835 


5.756 


11.28 


9.371 


8.402 


9.959 

9.219 


TABLE  VI 


Values  of  the  Parameter  i'CM^,  y) 
(Laminar  Nozzle  Heat  Transfer  Calculation) 


2 

3 

4 

5 

6 

1.2 

1.102 

.706 

.4505 

.3063 

.2235 

1.3 

1.069 

.7056 

.4705 

.3317 

.2473 

1.4 

1.04 

.701 

.4826 

.3498 

.266 

1.5 

1.012 

.694 

.4894 

.3626 

.2804 

1.6 

.9864 

.6856 

.4927 

.3713 

.2913 

1.667 

.9702 

.6795 



.4935 

.3755  | .297 

TABLE  VII 

Heating/Cooling  Efficiency  Q/<*  = S (Moo,  y) 
(Laminar  Nozzle  Heat  Transfer  Calculation) 
(assuming  <j  = 0.75) 




V 

2 

3 

4 

5 

6 

I 1.2 

10.58911 

10.47221 

7.742076 

5.201797 

3.407888 

1.3 

3.810589 

2.838522 

1.843574 

1.181419 

.7736683 

1.4 

2.069132 

1.353306 

.8382526 

.532963 

.3529183 

1 1.5 

1.345312 

.8203605 

.4981618 

.3174406 

.2127716 

1.6 

.9691441 

.565502 

.3403711 

.2178862 

.1476262 

{ 1.667 

.8067598 

.4611308 

.2768253 

.1778572 

.1212345 
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TABLE  VIII 


Values  of  the  Parameter  i"  (M,*,,  y) 
(Turbulent  Nozzle  Friction  Calculation) 


0.6617 


.6022 


.5508 


.506 


.4669 
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3 

4 

5 

6 

.324 

.1435 

.07301 

.04469 

.3078 

. 1474 

.0787 

.04824 

.2897 

.1465 

.08151 

.05072 

.2715 

.1429 

.08217 

.05205 

.254 

.1379 

.08138 

.05243 

.243 

.134 

.08031 



.05227 

TABLE  IX 

Values  of  the  Parameter  N (Mco,  y) 
(Turbulent  Nozzle  Friction  Calculation) 


TABLE  X 

Values  of  the  Parameter  G (Moo,  y) 
(Turbulent  Nozzle  Friction  Calculation) 


1.4 

1.5 

1.6 
1.667 


2 

3 

.06909 

.09623 

.06542 

.0836 

.0622 

.07464 

.05934 

.06779 

.05682 

.06238 

.05526 

.05933 

.08876 

~07561 

.06644 

.0616 


.1107  i 

.08717  j" 

.07259  j 

* — 4>- 

.06554 


.143 

.1028 

708069 

.07076' 


TABLE  XI 

Heating/Cooling  Efficiency  Q/^, 
(Turbulent  Nozzle  Heat  Transfer  Calculation) 


1.882194 


1.211807 


1.016553 


.9570714 


.9236017 


1. 

3 

1.42432 

.9746492 

.8367381 

.7888347 

.7648494 

1. 

4 

1.193997 

.8557685 

.749112 

.708328 

"76872611 

1. 

5 

1.056738 

.7847835 

.697265 

.6619956 

.6431093 

1. 

6 

.8855203  I 

i 

.738219 

.6628954 

.6317148 

.6148905 

1. 

667 

.9179515 

.7141839 

.6463789  ' 

.6170509 

.600834 

41 


1.667 


.6 


1.667 


TABLE  XII 


Values  of  the  Parameter  y) 

(Turbulent  Nozzle  Heat  Transfer  Calculation) 


.421 


.407 


.3924 


.3781 


.3344 


.3555 


.1709 


.1761 


.1778 


.1773 


.1754 

Tl736 


.07238 


.08028 


.08568 


.08912 


.0911 

J09181* 


ABLE  XIII 


alues  of  the  Parameter  K (M^,  Y,  Re^ 
bulent  Wake  Thickness  Calculation 


.5391567 


.5245009 


.5115746 


.4996897 


.4889489 


.4822006 


.6770355 


.6311099 


.5962623 


.5683407 


.5450785 


.5317594 


.8655276 


.7473666 


.6735038 


.6216873 


.5827305 


.5610909 


5 

6 

1.193497 

1.739235 

.9178681 

~T7l55874 

™ 762537“ 

7905835 

76887616  " 

.7680517 

.6285309 

'.6805207 

.5972383 

' T6372909' 

FIGURE  4.  The  function  S indicative  of  the  ease  with 
which  heat  is  exchanged  between  cusp  and 
laminar  nozzle  flow 
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FIGURE  5.  The  function  G indicative  of  the  total  cusp 
drag  in  turbulent  flow 
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COOLING/HEATING  "EFFICIENCY" 

OF  GDL  NOZZLES  IN  TURBULENT  FLOW 


POWER  EXCHANGED 


FIGURE  6.  The  function  Q /a  indicative  of  the  ease  with 
which  heat  is  exchanged  between  cusp  and 
turbulent  nozzle  flow 
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FIGURE  7a.  Laminar  lateral  profiles 
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DENSITY  DEFECT  FOR  Mm  = 4,  7 =1.4 
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FIGURE  10.  Density  defect  for  M®  = 4,  y - 1*4 


J 


wmmo'W  — 


TTTTTTTTTTTTTTrrrrnTTTTTTTTTTTTTTTTTTTTT 


TOTAL  TEMPERATURE  DEFECT  FOR  M = 4,  T = 1.4 
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FIGURE  11.  Total  temperature  defect  for  Ho,  = 4,  y 
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